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Abstract 

Fractal dimensions of eigenfunctions for various critical random matrix ensembles are investigated 
in perturbation series in the regimes of strong and weak multifractality. In both regimes we obtain 
expressions similar to those of the critical banded random matrix ensemble extensively discussed 
in the literature. For certain ensembles, the leading-order term for weak multifractality can be 
calculated within standard perturbation theory. For other models such a direct approach requires 
modifications which are briefly discussed. Our analytical formulas are in good agreement with 
numerical calculations. 
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I. INTRODUCTION 



One of the main results of quantum chaos is the statement that eigenvalues and eigen- 
functions of different physical quantum systems in the semiclassical limit are well described 
by universal ensembles of random matrices. When the corresponding classical model is in- 
tegrable, the Berry- Tabor conjecture [l| says that, once properly rescaled by the mean level 
density, quantum eigenvalues obey the Poisson statistics, i.e. they behave as eigenvalues of 
diagonal random matrices with independent diagonal elements. According to the Bohigas- 
Giannoni-Schmit conjecture quantum eigenenergies of classically chaotic systems on the 
other hand are distributed as eigenvalues of the standard ensembles of random matrix the- 
ory (RMT), where all matrix elements are independent random variables and the measure is 
invariant with respect to the conjugation over orthogonal, unitary, or symplectic matrices, 
depending only on the symmetries of the system (3j. A similar dichotomy exists between 
eigenfunctions. For classically integrable systems eigenf unctions are strongly localized on 
quantized tori. For classically chaotic systems eigenfunctions are fully extended and can be 
approximated by a superposition of elementary solutions with coefficients being independent 



random variables 



4]- 



Naturally, there exist systems which are neither chaotic nor integrable, and their random 
matrix description, if any, is not an easy problem. The tight-binding Anderson model 
in three dimensions with on-site disorder and nearest-neighbor coupling j^J is one of the 
most profoundly investigated examples. When disorder is not too large, its eigenfunctions 
are either localized or extended, depending on whether the corresponding eigenenergy is 
smaller or larger than an energy E c > called the mobility edge, which depends on the 
strength of the disorder. Eigenstates counted from the center of the band with \E\ > E c 
are exponentially localized and spectral statistics is close to the Poisson statistics. When 
\E\ < E c states are extended and the statistical properties of eigenvalues are well described 
by the usual RMT statistics. A new phenomenon appears when \E\ « E c , called the metal- 



insulator transition (MIT) point. In \&\ it was demonstrated numerically that statistical 
properties of eigenvalues and eigenfunctions around this critical point differ considerably 
from the standard statistics mentioned above. In particular, critical eigenfunctions are 
neither localized nor extended but have multifractal features, which manifest in a non-trivial 
scaling of mean eigenfunction moments (Y2i=i \ ^i\ 2q ) with the system size L. For localized 



states, moments do not depend on L, while for extended states they scale as L~( q ~ l > d , where d 
is the dimension of the system. For critical states the scaling is L~^ q ~ 1 ' Dq , where D q 6 [0, d], 
called the generalized (or multifractal) dimensions, are non-trivial functions of q. 

Numerical calculations for the Anderson model are limited by the fact that the system is 
three-dimensional. In order to investigate properties of critical states, simpler models have 
been proposed (see e.g. [7, 8]) in the form of N x iV random matrices with elements slowly 
decreasing away from the main diagonal, such as 

M mn = p m 5 mn + V(m - n) (1) 

where diagonal elements p m are independent random variables and off-diagonal elements 
V(m — n) decrease as the first power of the distance from the diagonal 

V(m-n) ~ - — 9 - — - (2) 
|m-n|>l \m — n\ 

The mostly investigated model of this type is the critical banded random matrix ensemble 
(CrBRME) of real symmetric (/3 = 1) or complex Hermitian (/3 = 2) matrices jsj, whose 
elements M mn are independently distributed Gaussian variables with zero mean and variance 
given by 

(|M nn | 2 ) ' 



P 

(\Mmn\ 2 ) = \ 



1 + 



— n \ 2 



m — n 
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m 7^ n. (3) 



The results obtained for this model are reviewed in [9|. For all values of the coupling constant 
g its eigenf unctions have nontrivial multifractal properties. Though fractal dimensions are 
not yet accessible to full analytical calculations, the construction of perturbation series yields 
an analytical approach to them. Perturbation series expansions of the D q were obtained for 
small and large values of the coupling constant g jslQ. For g 1 ("strong multifractality" 
limit, that is, the limit where states are almost localized) and q > 1/2, the leading term of 
the expansion into powers of g for eigenstates close to the center of the spectrum E = was 
found to be 

DqXi r±zm\T* " = 1 (4) 
^ T{q) I 575. " = 2 - 



When g 3> 1 ("weak multifractality" limit, where states are almost extended), the expansion 
of D q into inverse powers of g for eigenstates close to the center of the spectrum reads 

^ = 1 -ff^ + °fiV ( 5 ) 



2tt^ ' \g* 

Here we consider a different class of critical random matrix models, introduced in [12J. They 
are constructed from Lax matrices of classical integrable one-dimensional systems of N in- 
teracting particles whose positions and momenta are random variables. Integrability of the 
underlying systems makes it possible to analytically calculate various spectral properties, 
such as their joint distribution of eigenvalues and (in certain cases) several spectral correla- 
tion functions. A detailed investigation of spectral properties of these models was presented 



in 13]. The expressions obtained for spectral correlation functions have unusual features 
but their properties are characteristic of critical systems, therefore it might be expected that 
eigenf unctions display multifractality properties. The purpose of this paper is to investigate 
numerically fractal properties of several different critical random matrix ensembles, and to 
calculate analytically multifractal dimensions in perturbation series in the two regimes of 
strong and weak multifractality. Some of the results have been briefly mentioned in [fjj]. 

Let us introduce the four random matrix ensembles that we consider in the present 
paper. In all models, g is a free parameter (coupling constant) independent on N. The first 
three ensembles are related with Lax matrices of the rational, hyperbolic, and trigonometric 



Calogero-Moser models, respectively 



15]. The rational Calogero-Moser ensemble CM r is 



defined by N x N Hermitian matrices of the form 

M mn = p m S mn + ig- — — (6) 
m — n 

(S mn is the Kronecker delta). The hyperbolic Calogero-Moser ensemble CM^ consists of 
matrices 

with fi a real parameter independent on N. The trigonometric Calogero-Moser ensemble 
CM t is a set of matrices of the form 

Mr i • M(l fimn) / D \ 

mn = PmOmn + Iff — r , \~~7a7T ■ (8) 

JM sin [nym — nJ/A/J 

In all these ensembles, p m are independent random variables with zero mean and unit vari- 
ance 



CXXiAO. VV J 

HQ. 



The fourth ensemble that we consider is the Ruijsenaars-Schneider (RS) ensemble, defined 
as the ensemble of N x iV unitary matrices of the form 

with <3> m independent random phases uniformly distributed between and 2ir. It is related 
with the Ruijsenaars-Schneider model of N classical particles [l6] and with the quantum 
map corresponding to the quantization of an interval-exchange map 

Eigenfunctions and eigenvalues of matrices M mn are defined by the equation 

N 

M mn m n {a) = A a ^ m (a) (10) 

n=l 

(here and below we label eigenvalues and eigenfunctions by Greek letters). We assume that 
eigenfunctions are normalized, Yln=\ \^n{®)\ 2 = 1- The main object of our investigation is 
the asymptotic behavior of the mean moments of eigenfunctions 

1 N 

J ? = ^E P ^')' ( n ) 

i=i 

where 

1 N 

P*U) = T7^y( E \*M\ 2qs ( E - A -)) ( 12 ) 

is the local mean gth moment of eigenvalues. Here p(E) is the total mean eigenvalue density 

1 N 

= n(T, 6 ( e ~ A -)) ( 13 ) 

a=l 

and (. . .) denotes the average over the random matrix ensemble under consideration. Multi- 
fractal exponents D q characterize the asymptotic behavior of the mean moments when the 
matrix dimension iV goes to infinity. They are defined by the scaling 

I„ ~ N~ iq - l)D ". (14) 

For localized states the mean moments do not depend on N and thus D q = 0. For extended 
states D q equals the dimension of the system; in particular for matrix models eigenvectors 
are one-dimensional, and therefore D q = 1. Our purpose is to calculate the leading-order 
term of the perturbation series expansion of D q . We consider two regimes, characterized 
by strong and weak multifractality of the eigenfunctions. In the strong multifractality case 



the unperturbed states are localized and the zeroth order of the multifractal dimension is 
Dq°^ = 0. In the weak multifractality case the unperturbed states are extended and = 1. 
When a perturbation is added, fractal dimensions are changed to D q = D q ^ + d q . The small 
correction d q is obtained by expanding eigenfunction moments I q in perturbation series: 
when d q In N <C 1 , one has 

I q ~ N-( q -V D « } (1 - (q - l)d q \nN). (15) 

Perturbation series at small values of the coupling constant corresponds to strong mul- 
ractality and its construction follows the same scheme as for the CrBRME ensemble 



3, |9_|, [l0|, ll9[ . Matrices of CM ensembles can be expressed as 

M mn = p m 5 mn + g{\ - S mn )M^ n , (16) 

where diagonal elements p m are independent random variables distributed according to some 
probability density &{p). For matrices of the RS ensemble it is convenient to first rewrite 
the matrix (Q as 

e i$ m +j 7 r g (i-i/iV) sin^) e ^m 1-e 2 ^ 

M mn -d mn — s in(7T^/A^) + ^ ^ N 1 - &™{rn-n+g)fN ' ^ 

When g — Y the limits of each term exists and after rescaling by a factor exp[— brg(l — 1/N)] 
the matrix reduces at leading order to 

M mn * e^5 mn - ,(1 - S mn ) 2 fe^ i _ J [m _ n)/N , (18) 

which is of the form of (TIB"]) . The expansions ( IT6l) and ( IT8l) corresponds to the quasi- 
degenerateperturbation series, as different diagonal terms may be close to each other. Fol- 



lowing 
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19] , we obtain the first order of the multifractal dimension expansion into 



powers of g by replacing the N x N matrix M by all possible 2x2 submatrices of M 

(19) 



and summing over all indices m and n. The calculations in this regime are presented for our 
four models in Section [HI 

The situation in the regime of weak multifractality, where the unperturbed states are 
extended, is less simple. The CrBRME result ()5]) was derived by mapping the problem 



to the supermatrix sigma model. For our models (and other similar models as well) this 
approach seems not applicable and we use a more direct method. For the RS ensemble 
([H]), it is possible to construct a formal perturbation series expansion around any non-zero 
integer value of g. In Section II III we present detailed calculations of fractal dimensions for 
eigenvectors of this model at second order in e for g = 1 + e and in Appendix we briefly 
discuss the case g = k + e with integer k > 2. Fractal dimensions are an asymptotic property 
of the matrix ensemble when N — > oo, while perturbation series expansion requires to take 
e — > 0. It appears that in the case g = k + e with k > 2, contrary to the case k — 1, limits 
e — > and N — > oo do not commute and terms proportional to N appear in the expansion of 
eigenvalue moments (TO)]) , which are expected to be canceled by contributions of higher order 
in e. We postpone detailed discussion of this cancellation mechanism to another publication 



2CH ] . Here, in the case k > 2, we just extract the term proportional to IniV and neglect 
the contribution linear in N. The obtained results are in a good agreement with numerical 
simulations. 

For the Calogero-Moser ensembles (]6])-(j8]), the large g limit is obtained by dividing ( IT6|) 
by g and rescaling the variables, so that the models considered can be rewritten in the form 



M 



1 1 i-n 



M£l + -PmSmn, (20) 
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where Mml is a simple matrix with, in general, non-degenerate spectrum. For large g 
the second term can be considered as a perturbation and we use usual perturbation series 
formulas to calculate the fractal dimensions. The success of this approach depends to a large 
extent on the analytical accessibility of eigenvalues and eigenf unctions of the unperturbed 
matrix Mmli- In Section HVl we calculate the asymptotic expansion of the fractal dimensions 
for large values of the coupling constant for eigenvectors of CM r , CM/,, and CMj ensembles. 
It appears that for these ensembles, at second order in 1/g, fractal dimensions are zero for 
states close to the center of the spectrum. Numerical results suggest that for all three CM 
ensembles fractal dimensions decrease exponentially with g and therefore should be zero at 
any order of perturbation series in 1/g. By contrast, the perturbation series for states close 
to the spectral ends diverges for large matrix dimensions. Again, this corresponds to the fact 
that limits (7—7-00 and N — > 00 do not commute, and may be explained by the localization 



of such edge states at large coupling constant 



20] 



II. STRONG MULTIFRACTALITY 



19] and apply it to our systems 



In this section we recall the procedure used in 

Bay 

in order to derive the perturbation series expansion of the multifractal dimensions for small 
values of the coupling constant g. We present the method for the Calogero-Moser ensembles 
AS!)"©- Its adaptation to the RS ensemble is straightforward. 

Let us consider a matrix ensemble of the form ffl6l) . The first-order correction is obtained 
by considering only 2x2 submatrices of M in the sum PP. Setting h = gM^l, these 
submatrices are of the form 

Pm h 



(21) 



Eigenvalues of this matrix are given by 

M± = £± Vtf + W (22) 

with £ = (p m +p n )/2 and rj = {p m —p n )/2. The corresponding eigenvectors are (u + ,v + ) and 
(u-,V-) with 

h _ -A T 



and A-j. = rj ± \fff + \h\ 2 . Taking into account only the 2x2 submatrix (I2TT) in ( !T2|) means 
that the sum runs only over the two eigenvalues /i-t, so at leading order the correction to 
the contribution to P q {m) due to transitions to the state labelled by n is 

-L-^\ U+ \^5(E - p+) -6(E-t-r,) + \u_\ 2 «5(E - //_)}, (24) 

where the term 5 (E — £ — 77) corresponds to the zeroth order term. By assumption, h is 
fixed, and ensemble average is obtained by averaging over random values of p m and p n with 
a probability distribution a{p). The expression (1241) becomes 

^ / <W(£ + rj)„{t - v) ( m ^ )q 8{E - to) -8(E-t- V ) 
+ nJT^ E - (25) 



(\h\ 2 +^ 2 + y 

(the factor 2 comes from the Jacobian of the change of variables from {p m ,Pn) to (£,77)). 
The integral over £ is easily performed. Changing variables by setting rj — \h\ sinht in ff25l) . 
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one gets 

2 \ h \ /"j u f(y{E -\h\e- l )a{E -\h\e l 



p(E) 



J dfcoahf M \i + e -2t )q -a{E)a{E-\h\{e t -e- t )) 



+ <r(£+|%«)<r(£+|%-') . ( 2( . } 



;i + e 2i )9 

for the leading correction to P q (m). The distribution a coincides with the zeroth order of 
the level density p since unperturbed matrices are diagonal with eigenvalues pj. Thus at 
first order in h the functions of the form a(E ± \h\<p(t)) in (125]) can be replaced by p(E) and 
the integration over t can be extended over the whole axis. Equation (|26|) reduces to 



where we have used a continuation of the known integral 



2l| 



/o cosh 2fe (t) r(26) 1 ; 

The mean moment I q defined by ffTTl) is then obtained by summing this result over n ^ m, 
recalling that h = gMmn- Up to the first order in g the mean moment is therefore given by 



v^r(g-|) 



I q = l-2gp(E) r^j S (29) 
with 

1 N 

S =(nJ2\ M ^\)- (30) 

m,n=l 

The average is taken over remaining random variables entering Mmli (if any). For all critical 
systems S diverges logarithmically when iV — > oo, as 

S = 2slnJV + 0(l), (31) 



where s is a non-zero constant which depends on the model [7]. Identifying the coefficient 
of IniV in f ll5p and (I2"9j) one concludes that at first order in g multifractal dimensions of 
eigenfunctions with energy close to E are given by 



D q = 4gp(E)s ffa) ■ (32) 



This expression is valid only in the region q > 1/2. In order to find an analytic expression 
for q < 1/2 we rely on a symmetry that has been observed for the multifractal dimensions of 



many critical systems [22]] . Namely, the anomalous exponents defined by A q = (D q — l)(q— 1) 
are symmetric with respect to q — 1/2. In other words, multifractal dimensions for q < 1/2 
are related to those for q > 1/2 by A q = Ai_ 9 . This leads to the following analytic expression 
for q < 1/2, 

D, = l + ^^ fl ^A^i-lJ. (33) 

The explicit values of constant s for different models can be calculated as follows. For 
the rational Calogero-Moser ensemble CM r one has |Mmn| = l/\m — n\. The sum over m, n 
yields 

N N-l 

m,n — l k=l 

thus s — 1. A similar sum appears in the CrBRME model (j3j). For this ensemble one has 

* J" mnl , (35) 
V 2|m — n| 

where z mn is a Gaussian random variable with zero mean and unit variance which is real 
for (3 — 1 and complex for (3 = 2. Therefore (|^|) = \/2/ir for (3 = 1 and y / 7r /2 for (3 = 2, 



and summation over m, n yields, using ( 134"]) . s = l/y/n for /3 = 1 and s = a/vt/8 for = 2, 
which allows to recover (TjQ) with p(.E') — p(0) = 1/ a/27t. For the hyperbolic Calogero-Moser 
ensemble CM/j, one has 

1 ^ 9 JV ~ 1 AT I 

m,n=l fe=l ' 7 

and, as for CM r , s = 1. For the trigonometric Calogero-Moser ensemble CM t one has 

N „ 

.V 



1 ^ O N ~ l AT 1 

V I mnl AT 2 ^ |sin(/iA;/Ar)|' 1 ; 



m,n — 1 



The main contribution to the sum (I3"7j) comes from terms with k/N close to zeros of | sin(//x)|. 
The asymptotic behavior of the sum depends on the number of zeros when < x < 1. One 
obtains s = [fi/ir], where [.] denotes the integer part. For the RS ensemble Mmh is given 
by jUj) and it is easy to check that s = 1. 
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In order to assess the validity of the perturbation expansion formulas ( 132|) - (!33|) for strong 
multifractality, we performed detailed numerical calculations of multifractal dimensions for 
our four ensembles. Random realizations of matrices (ElHQ are diagonalized, and a fit of 
the moments of the form log(^ l^l 29 ) = a + blogN + c/N is obtained. For negative values 
of q, in order to avoid divergences due to exceedingly small values of the eigenf unction, 
a coarse-graining is first performed and the same fit as above is obtained for quantities 
l°s(Si(Sj For CM ensembles, the average is performed over eigenvectors around 

the eigenvalue E = 0, so that the density p(E) is assumed to be constant over the range 
of vectors considered. For RS ensemble the average is performed over all eigenvectors. 
In Figs. dHU we compare the numerical results for random matrix ensembles CM r , CM^, 
CM t and RS to formulas (13"2j) and (1551) . Without any fitting parameter, the agreement is 
remarkable. 



III. WEAK MULTIFRACTALITY FOR RS ENSEMBLE 

We now consider the opposite regime of weakly multifractal states. For RS ensemble 
this regime is reached when the coupling constant is g = k + e with e € 1 and k is any 
non-zero integer. Indeed, rewriting in this case matrix (Q as in ([1] 



e i$ m+i7 re(l-l/A0 sin ^ J*^ 1-e 2 ^ 

M mn - d m+ktn — s i n ( vre /iV) + ^ <Wm) n l _ e 2^{m-n+k+e)/N 

(here and below Kronecker symbols are to be understood modulo N), one sees that both 
terms have well-defined limits when e — > 0. The first term includes the constant phase factor 
e i7re(i-i/Af) Rather than taking into account the expansion of this term, it is more convenient 
to redefine matrix (jUJ) to M mn = M nn e~ me<yl ~ l / N \ Now one gets 

M mn = Mffi + eM« + 0(e 2 ) (39) 

where 

M^ = e^5 m+Kn (40) 

and 

(1) = _2vn i$m l-5 m+Kn 

lv± mn jy c I _ e 2-Ki{m+k-n)/N' v ' 

In this section we consider the case k = 1; the general case will be treated in the Appendix. 
For k — 1, eigenfunctions u m (a) and eigenvalues X a of the unperturbed matrix Mmn obey 
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FIG. 1. (Color online) Fractal dimensions D q as a function of q for CM r ensemble for g = 0.005 
(black circles), 0.025 (red squares) and 0.05 (green triangles). The pk are independent random 
variables distributed according to a Gaussian with mean and variance 1. Symbols are numerical 
results (symbols are larger than error bars), solid lines correspond to the formulas ()32p and (|33p 
with p(E) = \ j\p2/n and 8 = 1. Matrix sizes for numerical fit are N = 2 n , 8 < n < 13. Average is 
performed over the N/16 eigenvectors closest to the eigenvalue E = 0. Number of random realiza- 
tions of the matrix is between 2560 for N = 2 s and 40 for 2 13 . Dotted vertical line corresponds to 



Taking the product of both sides of this equation over all m we see that eigenvalues can be 
chosen as 



the equation 




(42) 



A, 



(43) 



a 
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FIG. 2. (Color online) Same as Fig. Q] for CM^ ensemble with \x = 2tt. The formula of Eqs. ([32 
and (I33j) is the same as for ensemble CM r (s = 1). 



where $ is the mean value of $j 



* = ^E%- ( 44 ) 

J'=0 

Eigenfunctions of the unperturbed matrix are obtained recursively from f|42|) by fixing one 
component, say u (a) = 1/yN, so that 

u n {a) = 7^ eiSnW ( 45 ) 



with 

S n {a) = - % + n$ . (46) 

i=o 

Contrary to the case A; = where unperturbed eigenfunctions are localized, for k > 1 
unperturbed eigenfunctions w n (a) are extended. 

When g = 1 + e, the perturbation of the matrix (Q at first order in e is given by fHTj) 

13 




FIG. 3. (Color online) Same as Fig. [T] for CMj ensemble with /i = 2ir. Here matrix sizes are 
N = 2 n + 1, 8 < n < 13. The formula of Eqs. ([32]) and flU is the same as for model CM r but with 
s = [ M /tt] = 2. 

with k — 1. Let us expand exact eigenfunctions in a series of unperturbed ones 

* n (a) = u„(a) + ^ C a pu n (f3). (47) 

Coefficients C a/ 3 can be expanded into series of e. The expansion of C a p with (3 ^ a starts 
with the first order in e and C aa with the second order. The modulus square of the eigenvalue 
components is 

\^n(a)\ 2 = \u n (a)\ 2 + \<(a)C a pU n (/3) + u n (a)C* a(S u* n (P) 



where 



J2C a pC* aj u n (f3)u*M = 1(1 + A n (a) + B n (a)) 

0,-y 



A n (a) = [e iSnW - iSn{a) C a p + e -iSnW+iS n{ *) c ^ 

P 



(48) 



(49) 
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FIG. 4. (Color online) Fractal dimensions D q as function of q for RS ensemble for g = 0.01 (black 
circles), 0.05 (red squares), 0.1 (green triangles) and 0.2 (blue diamonds). The are independent 
random variables distributed uniformly in [0, 2ir]. Symbols are numerical results (symbols are larger 
than error bars), solid lines correspond to the formula of Eqs. (|32h and (I33f) with p(E) = 1/(2tt) 
and s = 1. Matrix sizes for numerical fit are iV = 2 n ,8<n<12. Average is performed over 
all eigenvectors. Number of random realizations of the matrix is from 32 for N = 2 8 to 2 for 2 12 . 
Dotted vertical line corresponds to q = \. 



and 




(50) 



0,7 

Taking (jlSjl to the gth power we get, up to the second order in e, 



*n(a)| 29 = ^ 1 + q(A n (a) + B n {a)) + S^Z 




(51) 



Normalization of the wavefunction implies that 




(52) 



n=l 



n=l 



15 



and therefore from (151]) we get 



N 



n=l 



N 



2N 



a 



n=l 



At leading order in e, coefficients C a p with /3 7^ a are given by 



C a p — e 



X a — A 



with 



Using (P]) - (H6]) we obtain 



= X)<(/3)M««n(a). 



A _ A = e i# +27 ria/iV _ e i*+2^/iV = % ^ * (q _ p^^a+fi/N 



and 



a/3 



— T 



jS n {a)-iS m {P)+i<b m -m{m+l-n)/N 



iV 2 ^ sinfvrfm + 1 - n)/j\fl 

n,m L v 



Consequently, 



me 



E 



M+iXmn (*) 



-7ri(m+l— ra)/iV 



where 



2iV 2 sin[7r(a; — /3)/N] j—^ sm[ic(m + I — n)/N] 



r, / \ 2tt . 1, 2tt„, 1, 
F Q/3 (m, n) = j^a{n - -) - + -) 



and the function x m „($) contains the dependence on random variables 

n— 1 m 

According to ( l53l) . the fractal dimensions are obtained by calculating 



+ ^$ i + ( n - m -i)$ . 

j=0 j=0 



1 * 

n=l /3,7 



where 



N 



n=l 



a/3 



v rp iS 'n(7)-i'S , „(a)^-y . -iS n (7)+iS„(a) f-r* 

a e <^a7 ~r e 
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Using expression (146!) for S n (a) and summing over n in (162]) one gets 

i?(a, /3, 7) = C aP (C* ai 5^ + C on 62 a -p,* f ) + c.c. (63) 

where c.c. denotes the complex conjugate. Finally at leading order in e 

W(a) = (\C^\ 2 + C a pC* a>2a _p +pN ) + c.c. (64) 

where p = or ±1 accounts for the fact that the delta function in (16"3"1) is modulo N. 

In order to obtain fractal dimensions averaged over all eigenvectors of our matrix we need 
to calculate 

a 

From ( )58|) we have 



lCa ^ 2 ~ 4N^m 2 [n(a-(3)/N] (66) 

e -iF a(3 (m',n')+i-F Q/3 (m,n)-ix m /„/(*)+iXmn(*) 



E 



^—m(m—m'—n+n')/N 



sin[7r(m / + 1 — n')/N\ sin[7r(m + 1 — n)/N\ 

n' ,m' ,n,m 
n 1 y^m' +l,n^m+l 

The sum of |C a/ 3| 2 over a and (3 involves terms of the form 

-iF a f,(m' ,n')+iF ol 8(m,n) N ~ 1 „27ri/3 (m'-m)/iV 

e P P _ ST^ e 27ria(n-n'+m'-m)/N s^s 

sin 2 [vr(a- 0)/N] ^ sin 2 [7r/3/An ^ 1 ' 

which are non-zero if and only if m — n = m' — n' . We then take the average of (1661) over 
random phases uniformly distributed in [0, 27r]. One has to consider averaged quantities 
of the form 

which are non-zero if and only if coefficients of all the $j in the exponential are zero. The 
condition m — n = m' — n' cancels out terms in $ in ( |68|) , and we have 

~ Xm>n>(*) + Xmn(*) = ^ / j , " (69) 

Note that the case m = n — 1 is excluded since it does not appear in the sum ( )66|) . In 
all cases, the only way to have all $j vanish from (169]) is to have m = m' and n = n! . 
Performing the remaining sum in (|66l) under this condition, we get 

( ^ ? |Cq/?|2) = 4A^ ? sin 2 (7r/3/iV) ^ sin 2 (vrn/iV) ' (?0) 

« /3=1 v ' / / n= l \ i i 
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The second term in (1M1) can be calculated in the same way, using 

_ f -7ri(m+l-n)/AT 

o Q ,2a-^+pjv - m2 sin ^ a _ ^ ^ s in[7r(m + 1 - n)/N] 

(71) 

The sum of C*^C*^ 2a- ( g+pAr over a and /3 involves terms of the form 

ST^ ST^ e P 2 2 = e 27ria(n+n'-m'- m -2)/yV 

sin 2 [vr(a-/3)/iV] " sin 2 [vr/3/iV] ^ 

(72) 

which are non-zero if and only if n + n' = ml + m + 2. Again, when averaging over random 
phases, this latter condition cancels out terms in $, and 

- w.) - x m „(*) = < - 5" *' + 5r 1 * J m - " m 



All $j terms vanish in the expressions above if and only if m! = n — 1 and n' = m + 1. 
Finally 



r 2 £ 2 ^V- 1 j N - 1 e -2iTv/3n/N 



(E E C ^ C ^a-fi+^) ~ 4A T2 E sin 2 (7r/ 3 /iV) E sin 2 (7m/i V) • (74) 

a /^a /3=1 \ ' / / n= l \ / / 

Inserting ( J70|) and ( 1741) into f lMl) and performing the sum over /3 using the following identity 
valid for integer n, 

E 1 sin 2 (vT/3n/iV) 
A- K / = n(N-n), (75) 
^ sin 2 (7r/3/iV) 1 J ' 1 J 



we obtain the final expression 



2 2 N-l 



(W) = —Y <N ~ n) (76) 
^ 3 ^sin 2 (7rn/A^) 1 ' 

The summand in this expression diverges when N — > oo when n is close to or N. One can 
rewrite the sum as 

n^ n(N~n) _n^( n(N - n) N N \ ™ 1 

N 3 ^sm 2 (7rn/N) N ^ \ N 2 sin 2 (vm/iV) vr 2 n vr 2 (iV-nW ^n' 1 J 

n=l v ' ' n=l x \ i / \ i / n=l 

For N — > oo, the first sum tends to a finite value 



VT 2 



f ( ^t4 - 4 277^) d y = 2[1 - M2tt)], (78) 

Jo Vsm 2 (vT?/) 7T 2 ?/ vr 2 (l-?/)y 
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while the second sum diverges as In N. Thus when N — y oo 



(W) = 2e 2 lnA + O(l). (79) 

From ([S3]) and (T6Tj) it follows that 

jyBE i*»i 29 > ^ iVl ^( 1 + ?(9 - 1 > 2 lniV ) ( 8 °) 



a n=l 

from which one extracts the multifractal exponents D q using (I15j) (recall that g = 1 + e) 

D 9 = 1 - q(g - l) 2 . (81) 

In the Appendix this result is generalized to g = k + e with integer k > 2. It is argued that 
fractal dimensions at leading order are then given by 

D^i-^ 1 ^ 1 - (82) 

In Fig. [5] we present the result of numerical calculations for D\ in the RS model for different 
g. The agreement with the above perturbation series results is quite good. Numerical 
results suggest that flHTl) is actually the exact formula for the multifractal dimension D\ of 
the RS ensemble in the whole range < g < 1. The calculation of higher-order terms of 
the perturbation series, though in principle possible, may be more complicated. The point 
is that in the limit iV — y oo and fixed e spectral properties of the RS ensemble are different 



for e > and e < 



12 



13]. If the conjecture in [14j is correct, it will signify that fractal 
dimensions for this ensemble also depend on the sign of e, and simple analytical results ( 18 ip 
and (|82|) may be valid only at second order in e. 



IV. WEAK MULTIFRACTALITY FOR CM MODELS 

We now investigate fractal dimensions for the CM ensembles when the coupling constant 
g goes to infinity. Matrices of these ensembles are defined as the sum of a diagonal term 
and a Hermitian off-diagonal matrix (see (E])— (JH])) • In order to apply the perturbation series 
analysis for large g we redefine our matrices as M mn = M mn /g, so that 

+ -PmSmn, (83) 
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FIG. 5. (Color online) D\ as function of g for the RS model. Circles are numerical results obtained 
by averaging over all eigenvectors taken from 128 realizations for N = 2 8 to 8 realizations for 
N = 2 12 . Dashed lines indicate the leading perturbation series results (|82p for q = 1. Solid line is 



an eye guide corresponding to the theoretical value given in [13 ] . assuming that the conjecture of 
\l\ holds. Inset: same data zoomed in for 1 < g < 3. 



where Mi°| are the off-diagonal terms in ©-([E]) (divided by g). These matrices are non- 
random Toeplitz matrices, as their entries only depend on the difference m — n. The prop- 
erties of such matrices when iV — > oo are well investigated (see e.g. {23] and references 
therein) . 

Exact eigenvalues and eigenvectors of a Toeplitz matrix are not known in general. In the 
case of the CM t ensemble with /i = 2-k the matrix Mmn becomes a circulant matrix, and 
as such its eigenvectors are simply given by Uk(ct) = exp(2'nrak / N) / \/N . The situation is 
thus similar to that of the previous Section (cf. (H5|) ). where unperturbed eigenvectors are 
extended. We now consider this case by closely following the approach of the RS ensemble. 
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The unperturbed matrix reads 

M(°) = 2nKl - 6mn) . (84) 

mn N sin [2vr {m-n)/N\ K ' 

For simplicity we choose N = 2K+1 as an odd integer. It is convenient to define eigenvectors 
in the following manner 

u n (a) = -JL exp (j^n a) (85) 

with index a running from — K to K. With such a choice, the corresponding eigenvalues 
take the simple form 

9 JV ~ 1 

Aa " " iV ^ sin(27rn/iV) e " " iV " (86) 

n=l v ' ' 

The spectrum is rigid, which implies that the density of unperturbed eigenvalues is constant 
in the interval (— 27r,27r). All equations from ( H71) to ( 155]) remain valid, with e = 1/g, 
Mmli = Pm$mn and S n (a) = Airna/N. From Eq. (JSJ) one obtains 

1 N 

r a - T> p 47rim ( a -^/ N rs7^ 

Equations (l6Tj) -f l64|) are unchanged, thus we just have to calculate W(a) from (jSlj) and 
(!H7|) . In both terms appearing in f|64l) . the averaging over random variables p m involves 
expressions of the form 

N 

{PmPm>}e ±4ni{m - m ' ){a -P )/N = N{p 2 }. (88) 

m,m'=l 



Thus after averaging over p m (IMl) gives 

W«» = S E 2A ° - ^ - ^i£? . (89) 

Since indices a and /3 run between —K and i^, we have 2a — (3 £ [— 3i^, 3iT]. Replacing the 
unperturbed eigenvalues by their expression ( |86l) we get 2A a — \p — \2a~/3+ P N = ^iip, which 
means that for 2a — (5 £ [— K, K] (i.e. p = 0) the contribution to W vanishes. Consequently, 
contributions to ( |64l) come only from terms with > 2a + K or (3 < 2a — K, yielding 



N 2 (p 2 ) 

(W{a)) = -^gj£T{a,N) (90) 
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where 

j=K+l-\a\ J V Jl 

This gives us the final expression for the averaged momenta of eigenvectors 

n=l ^ ^ 

While the RS ensemble is an ensemble of unitary matrices, matrices corresponding to CM 
ensembles are Hermitian, and the position of the eigenvalue in the spectrum, specified by 
the value of the index a G [— K, K], does matter. For eigenvectors whose eigenvalues are 
close to the center of the spectrum, the index is such that |a| <C K. For \a\ fixed one has 

T(a,JV)— >0, (93) 

N— >oo 

which, using ffl5|) . implies that second-order correction to the multifractal dimension is for 
these states. 

Numerical calculations, as presented in Fig. |6l suggest that in fact multifractal dimensions 
for Calogero-Moser models at large g are exponentially close to 1, for example, 

Di ~ 1 - C{g)e- b9 (94) 

with a certain slow varying function C(g) and some parameter b. If this is indeed the case, 
any order of perturbation series in 1/g would yield zero for states with \a\ <C K. 

By contrast, eigenfunctions corresponding to eigenvalues at the edge of the spectrum, 
labelled by a with |a| ~ K, are localized. Indeed, let v — K — \cx\ be the index of the 
eigenvector counted from the edge of the spectrum. One can show that for v fixed TiK — 
v, N) can be lower and upper bounded by functions which asymptotically behave as N/ 
and N/ v, respectively. The perturbation series expansion ( 192]) is valid formally only when the 
second term is much smaller than the first one, that is, only if iV £ where £ ~ g 2 u. Thus at 
fixed values of g and v ( 192]) is no longer a correct approximation for large N and calculation 
of higher-order terms becomes important. In other words, limits e — > and iV — > oo do 
not commute. The breakdown of the above perturbation series seems to indicate that our 
assumption that eigenfunctions with fixed u are almost extended is not correct. Numerical 
simulations show that for large g eigenvectors associated with eigenvalues located at the end 
of the spectrum are localized and that the localization length increases with the distance from 
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FIG. 6. (Color online) D\ as a function of g for models CMj (red circles), CM r (blue squares) 
and CM/, (black triangles), obtained from 256 realizations of the random matrices for N = 2 8 
to 32 for N = 2 11 (only N/16 central vectors are considered). Dashed lines indicate the leading 
perturbation series results for small g: D\ = 2\phisg with s = 1 for models CM r and CM/,, and 
s = 2 for model CMj. Solid lines indicate the fits for large g: D\ ~ 1 — 0.85e~ 5 ' 89 for model CM r , 
0.78e" 5 - 5 f for model CM h , and I>i 1 - .92e~ 12 - 2 f for model CM t . 



the spectral boundaries. Rough estimation gives that the localization length is proportional 
to £ = g 2 v. When ( C iV states are strongly localized and fractal dimensions are zero. 
But when £ exceeds N localization becomes unessential and fractal dimensions take non- 
trivial values as in (|94|) . This situation is analogous to the 3-dimensional Anderson model 
where states close to spectral boundaries are localized. The existence or absence of a sharp 
mobility edge for the Calogero-Moser models requires further investigations. A more careful 
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discussion of such boundary states in these and other models will be given elsewhere 

As an example, we display in Fig. [7] the inverse participation ratio, which is the usual 
measure of the localization length, averaged over eigenvectors taken at different positions in 
the spectrum. It is clearly seen that states close to spectral boundaries are strongly localized 
while states far from boundaries are delocalized. 

For models CM r and CM/,, where the unperturbed matrices are no longer circulant, there 
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FIG. 7. (Color online) Inverse participation ratio £ = l^i( a )| 4 ) as a function of matrix 

size N for model CMf for g = 0.5 and [i = 27r. The average is taken over 32 realizations of the 
random matrix for N = 2 8 , to 2 realizations for TV = 2 12 , and over eigenvectors associated with 
eigenvalues A2 r fc+i < ••• < A2 r (fc+i) with r = 8 and k = 0,1,2 from bottom to top (eigenvalues 
are ordered so that Ai < ... < Aj\r)- We choose p m as independent random variables distributed 
according to a Gaussian with mean and variance 1. Inset: same as the main panel (three lower 
curves) and eigenvectors averaged over Ajv/2-r> • • • ; ^N/2+r (black circles), Ajv/4_ r > • • • 3 AAr/4 +r (red 
squares), Ajy7 8 _ r; . . . , A^v/ 8+r (green triangles), r = 8. For the three upper curves, solid lines are 
linear fits with slopes respectively given by 0.99, 0.97, and 0.94. Logarithm is decimal. 

is (to our knowledge) no exact analytical expressions for the unperturbed eigenvectors and 
eigenvalues. Approximating eigenfunctions for these models as Fourier harmonics similar to 
(|85p and calculating eigenvalues as in (I5rj|) one finds a good agreement with the numerically 
computed spectrum. Using such an approximation we find that the behavior for large g in 
all Calogero-Moser models is very similar. Numerical results presented at Fig. [6] agree with 
this conclusion. 
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V. CONCLUSION 



We investigated fractal dimensions of eigenfunctions for critical random matrix ensembles 
dEJ)-© m perturbation series for strong and weak multifractality. Spectral statistics of these 



ensembles are very different 12] but fractal dimensions have many common points. 

For strong multifractality, when the coupling constant is small, fractal dimensions have 
the universal form 



D q = 4gp(E) s 



r(?) 



(95) 



where p(E) is the density of diagonal elements, and s is a constant which depends on the 
system. For weak multifractality, fractal dimensions D q to leading order also have the 
universal form 

D q = l- qt, (96) 

with a constant t depending on the system. For the RS ensemble (jHJ) the weak multifractality 
regime corresponds to the vicinity of all integer points g = k with \k\ > 1, while for the 
CM ensembles it corresponds to g — > oo. The values of constants s and t for the models 
considered are presented in Table [B 



Model 


CrBR 

= 1 


ME 

P = 2 


RS 




CM r 


CM h 


s 






1 


TV 


1 


1 


t 


1/(2tt0) 


1/(4tt 5 ) 


(g-kf/k 2 












TABLE I. Constants s and t for different models. 



For the CM ensembles the second order contribution to fractal dimensions for states far 
from spectral boundaries is zero. Numerical calculations suggest that fractal dimensions for 
these ensembles are exponentially small for large g and, therefore, are zero in any finite order 
of perturbation series in 1/g. 

The methods used in the paper are general and can, in principle, be applied to a large 
variety of critical ensembles. Properties of strong multifractality are well understood and 
controlled but weak multifractality limit for certain ensembles requires further investigation 
20]. 
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Appendix A: Calculation of fractal dimensions in perturbation series for the RS 
ensemble in the general case 



The purpose of this appendix is to calculate multifractal dimensions for the RS ensemble 
to leading order in e when g = k + e for integer k > 2. The unperturbed matrix in this case 
is a matrix with shift by k, 

Afl = e** m <W*,n • (Al) 
Eigenfunctions u m (a) and eigenvalues \ a obey the equation 

e^ m u m+k (a) = X a u m (a) . (A2) 

If one component u m (a) is fixed, then Eq. (IA2j) allows to deduce all components of the form 
u m+t k((y). Thus, when gcd(k,N) is different from 1 one can construct families of solutions 
of the eigenvalue equation flA2[) whose coefficients are nonzero only on subsets of indices of 
the form t + rk, r = 0, 1,2, . . .. In what follows we denote c =gcd(k,N) and R = N/c. 
Indices n with < n < N — 1 are replaced by a pair of indices (t, r) with < t < c — 1 
and < r < R — 1 such that n = t + rk. Fixing u t fi(a) = l/\/R, the other components are 
obtained recursively from (1A2j) as 

u t , r (a) = -±=e iSt ^ (A3) 



with 

S t >) = -^r - ]T $ tJ + ri> t . (A4) 
i=o 



1 s^R-l 

values can be chosen as 



where $j is the mean value of $tj, defined by $4 = S^Lq 1 ®t,r ^ ne corresponding eigen- 



A t)Q = e i4,t+2 ™/ fl . (A5) 

The treatment is very similar to the case k = 1, the only difference being that families 
of eigenvectors with fixed t are treated in parallel. To simplify the discussion, let us now 
consider the case gcd(/c, N) = 1. The perturbation series of the matrix (Q up to the first 
order in e is given by fl39|) - fj4T|) . Equations P7|) - fl55|) are still valid and they yield 

lire J e iF Q/3 (m,n)+iXmn(*) 
C a = — — - - m (rn+l-n)k/N /A fi N 

ap 2A 2 sin[7r(a-/3)/A] ^ sin[vr(m + 1 - n)A;/A] ' 1 ' 
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where F a p(m,n) is given by (1591) . and 

n— 1 m 

Xm „($) = - ^ $ fcJ + ®k,3 + fa - m - 1)<£ . (A7) 
i=o j=o 

A straightforward generalization of all calculations done for the /c = 1 case gives the final 
expression for k > 1 

w < A8 > 

Terms divergent when N — > oo correspond to regions of n close to zeros of the denominator, 
n = mM + r where M = [iV/fc], m = 0, 1, . . . k, and |r| <C M. Terms with m = and m — k 
can be treated as in Section HTT1 and give a logarithmically divergent contribution 2ln(N)/k 2 
but contributions of regions with m = 1, . . . , k — 1 grow linearly with N. Their calculation 
can be done by extending the summation over r to all integers. Finally we gets 

(W) =aN + ^- In N + 0(1) (A9) 
k 2 

where 

7r 2 x ^ m(k — m) 

a = 7J> —T} 77T A10 

A; 4 ^ sin 2 Urns k) 

and s is the residue of A^ modulo k (we consider here the case of co-prime A^ and k). Taking 
into account the logarithmic term results in f|82|) . 

The existence of linear in A^ term, as above, indicates a breakdown of the simple pertur- 
bation series approach and the necessity of a more careful treatment of higher-order terms. 
The precise investigation of the compensation of such terms is beyond the scope of this 
paper and will be analyzed elsewhere j^. 
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